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Results are reported from a search for active to sterile neutrino oscillations in the MINOS long- 
baseline experiment, based on the observation of neutral-current neutrino interactions, from an 
exposure to the NuMI neutrino beam of 7.07 x 10 20 protons on target. A total of 802 neutral- 
current event candidates is observed in the Far Detector, compared to an expected number of 
754±28(stat) ±37(syst) for oscillations among three active flavors. The fraction f s of disappearing 
that may transition to v s is found to be less than 22% at the 90% C.L. 
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The disappearance of muon neutrinos as they propa- 
gate from their production source is well established by 
experimental evidence accumulated over the past several 
decades [l|-[5| . Although the deficit of charged-current 
(CC) interactions is generally interpreted as due to os- 
cillations between the weak flavor states of active neu- 
trinos, with — > v T transitions representing the dom- 
inant channel [4 6] , more exotic scenarios where active 
neutrinos oscillate into an unseen sterile neutrino flavor, 
u s , are not excluded. The possible existence of one or 
more light sterile neutrinos, in addition to the three ac- 
tive flavors, has been widely discussed Q and could con- 
tribute to the understanding of the neutrino mass spec- 
trum Q or explain apparent differences in behavior be- 
tween neutrino and antineutrino oscillations [j| . Interest 
in sterile neutrinos has been renewed with the latest ob- 
servations from antineutrino running in the MiniBooNE 
experiment, which may be explained by mixing models 
incorporating one or more sterile neutrinos ■ Further- 
more, recent results from the WMAP experiment might 
suggest the existence of a fourth neutrino generation with 
mass less than 0.58 eV 



MINOS can probe active to sterile neutrino mixing 
driven by the atmospheric mass-squared splitting by 
measuring the rate of neutral-current (NC) events at two 
locations, over a baseline of 735 km. Because NC cross- 
sections are identical among the three active flavors, NC 
event rates are unaffected by standard neutrino mixing. 
However, oscillations into a sterile nonintcracting neu- 
trino flavor would result in an energy-dependent deple- 
tion of NC events at the far site. This letter reports 
results from a search for sterile neutrino mixing, using a 
data sam ple twice as large as that used in previous pub- 
lications [l2j. The data are compared to models where 
these neutrino oscillations are driven either by the at- 
mospheric mass scale Am| 2 alone or along with a mass- 
squared splitting Am 2 3 having magnitude 0(1 eV 2 ). 

MINOS measures neutrinos from the NuMI beam [l3| 
using two detectors: the 980 ton (27 ton fiducial) Near 
Detector (ND), located 1.04km downstream of the beam 
target at Fermilab; and the 5.4kton (4.0kton fiducial) 
Far Detector (FD), placed 735 km downstream of the 
target in the Soudan Underground Laboratory, in Min- 
nesota [3] ■ The detectors are planar steel and scintillator 
tracking calorimeters. Each plane is composed of 2.54 cm 
thick steel and 1 cm thick plastic scintillator arranged 
in 4.1cm wide strips. The energy resolution function 
for neutrino-induced hadronic showers is approximately 
56% /VE [15J. The data were collected in an exposure 
of 7.07 x 10 20 protons on target taken exclusively with a 
beam configuration for which the peak neutrino event en- 
ergy is 3.3 GeV. The NuMI beam includes a 1.3% {v e +u e ) 
contamination primarily from the decay of muons origi- 
nating in kaon and pion decays. 

In addition to the increased statistics, the new analysis 
includes changes to the shower reconstruction, whereby 



active strips with fewer than two photoelectrons of pulse 
height are excluded from the clustering algorithms, re- 
ducing effects from crosstalk. The analysis also benefits 
from a complete reevaluation of both the event selection 
criteria and the effects of systematic uncertainties. 

The ND registers a high event rate during operation, 
with multiple neutrino interactions occurring through- 
out the detector for each beam spill. The total activity 
recorded during a spill is separated into activity slices us- 
ing timing and spatial criteria Ideally, each activity 
slice would correspond to one neutrino interaction, but 
some failure modes result in one activity slice containing 
information from different interactions, for which sepa- 
rate events may then be reconstructed. Simulations show 
that these failure modes increase the number of interac- 
tions selected as NC with reconstructed energy (E leco ) 
lower than 1 GeV by 37% in the ND. This background 
is reduced to 11% by removing a reconstructed event if 
it contains less than half of the total energy deposited 
in the activity slice, or if the event has fewer than three 
contiguous planes with at least two photoelectrons read 
out in each plane [HI]. 

Only a few beam-related events are recorded each day 
in the FD fiducial volume. Interactions are selected for 
the analysis if they occur between 2/xs before and 12 /is 
after the expected start time of the 10 /j,s spill at the FD. 
Possible backgrounds due to detector noise or cosmic-ray 
muons in coincidence with the spill window are removed 
by various selections 12| . The remaining nonbeam back- 
grounds after application of these criteria represent only 
0.5% of the expected NC interaction rate in the FD. 

In the MINOS detectors, NC interactions give rise to 
events with a short diffuse hadronic shower and either 
small or no tracks, whereas CC events typically display 
a long muon track accompanied by hadronic activity at 
the event vertex. Events crossing fewer than 47 planes 
for which no track is reconstructed are selected as NC; 
events crossing fewer than 47 planes that contain a track 
are classified as NC only if the track extends less than 
6 planes beyond the shower [16j | ■ These selections result 
in an NC-selected sample with 89% efficiency and 61% 
purity. Events that fail both NC criteria are selected as 
CC if they pass the classification procedures used by the 
MINOS muon neutrino CC disappearance analysis [H; 
otherwise they are removed from the analysis. Highly 
inelastic and CC events, where the muon track 
is not distinguishable from the hadronic shower, are the 
main source of background for the NC-selected spectrum. 

The predicted NC energy spectrum in the FD is ob- 
tained using the ND data. An estimate of the ratio of 
events in the FD and ND as a function of reconstructed 
energy, E rcco , is calculated from Monte Carlo simula- 
tions. The ratio is multiplied by the observed ND energy 
spectrum to produce the predicted FD spectrum [171 ]. 
Alternative methods of predicting the FD energy spec- 
trum 18, li| yielded very similar results. To avoid biases, 
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FIG. 1: Reconstructed energy spectrum of NC-selected events 
in the ND (solid points) compared to the Monte Carlo pre- 
diction (open histogram) shown with 1 a systematic errors 
(shaded band). Also displayed is the simulation of the back- 
ground from misidentified CC events (hatched histogram). 



FIG. 2: Reconstructed energy spectrum of NC-selected events 
in the FD (points with statistical errors) compared with pre- 
dictions for standard three-flavor mixing with and without 
u e appearance at the MINOS 90% C.L. limit [13] (dashed 
and solid lines respectively). 



the analysis selections and procedures were determined 
prior to examining the FD data. Figures [T] and [5] show 
the reconstructed energy spectra in each detector. The 
relatively low number of events observed in the ND en- 
ergy spectrum for E ICCO < 1 GeV is a consequence of the 
application of the ND-specific selections described above. 

The analysis classifies 97% of ^-induced CC events 
as NC, requiring the possibility of v e appearance to be 
considered when extracting results. The normal neutrino 
mass hierarchy is assumed with 0j3 = 11-5° and 5cp = ft 
at the MINOS 90% C.L. limit Efl. Mikheycv-Smirnov- 
Wolfcnstcin-like matter effects [21| , due to differences be- 
tween the matter potentials for active and sterile neutri- 
nos, are at the subpercent level for the MINOS baseline 
and are neglected in this analysis. 

The selection procedures identify 802 NC interaction 
candidates in the FD, with 754±28(stat)±37(syst) events 
expected from standard three-flavor mixing (assuming 
013 = 0°). An excess relative to the 6*13 = 0° predic- 
tion is observed in Fig. [2] for 1 < E rcco < 5 GeV. This 
excess does not present significant evidence for new neu- 
trino phenomena and is treated as a statistical fluctua- 
tion. The agreement between the observed and predicted 
NC spectra is quantified using the statistic R: 

n _ Vdata - B C C ,-, x 

H = o > i 1 ) 

where Adata is the observed number of events, Bqc is 
the predicted CC background from all flavors, and Snc 
is the expected number of NC interactions. The values 
of Adata, 5*nc and contributions to Bqc for various re- 
constructed energy ranges are shown in Table HI The 
values of R obtained for each energy range show no evi- 
dence of a depletion in the NC flux at the FD, support- 
ing the hypothesis that standard three-flavor oscillations 



explain the data. A value of R = 1.09 ± 0.06(stat) ± 
0.05(syst) — 0.08(v e ) is measured over the full energy 
range 0—120 GeV, where the last term is the change re- 
sulting from inclusion of v e appearance at its maximally- 
allowed value from the MINOS 90% C.L. limit. There- 
fore, the depletion of the total NC event rate is less than 
3.2% (11.2%) at 90% C.L., where the value in parentheses 
is obtained assuming v e appearance. 
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TABLE I: The R statistic and its constituent components 
for several reconstructed energy ranges. The numbers shown 
in parentheses include v e appearance with 613 = 11.5° and 
Sep = 7T. The displayed uncertainties are statistical, system- 
atic, and the uncertainty associated with v e appearance. 



The data are compared with two models of neutrino os- 
cillations that allow admixture with one sterile neutrino. 
In the first model, identified as = mi, the first and 
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fourth mass eigenstatcs are treated as degenerate and 
the oscillatory behavior is assumed to be driven only by 
the atmospheric mass scale. The second model, referred 
to as m 4 3> m 3 , assumes a large difference between the 
fourth and third mass eigenstates, introducing an addi- 
tional mass scale Am| 3 with magnitude 0(1 eV 2 ), so that 
no oscillation-induced change of the neutrino event rate 
is measurable at the ND site, but rapid oscillations are 
predicted at the FD location. The latter model is also 
sensitive to portions of the region of interest studied in 
the MiniBooNE antineutrino data [22| . Detailed descrip- 
tions of these models are provided in Refs. 12J and [23| . 

Both the NC-selected energy spectrum shown in Fig. [2] 
and the CC-selected spectrum in the FD data are used 
in the fits to the oscillation models. The dominant sys- 
tematic uncertainties, discussed below, are added as nui- 
sance parameters to the \ 2 function used in the fits. The 
best-fit values are summarized in Table UH which dis- 
plays results with null and maximal v e appearance. The 
best=fit value obtained for |Am| 2 | is consistent with the 
results from the muon neutrino CC disappearance anal- 
ysis Q. A 90% C.L. limit on the sterile mixing angle 
of 634 < 26° (37°) is found for the 7774 = mi model. 
For the 7774 ^ 7773 model, limits of 624 < 7° (8°) and 
6*34 < 26° (37°) are obtained at the 90% C.L. j2J|. Pure 

— > v a oscillations are excluded at 99.8% (96.2%) C.L. 
The numbers in parentheses represent the limits ex- 
tracted for maximal v e appearance. 

The coupling between active and sterile neutrinos may 
also be quantified in terms of the fraction of disappear- 
ing that oscillate into v s , f s = P V ^ V J(1 - P^^J, 
where the P Vli -*v x refer to neutrino oscillation probabil- 
ities. For both sterile oscillation models, the limit set 
on the disappearance fraction is f 3 < 0.22 (0.40) at the 
90% C.L. This limit represents a 57% improvement on 
the previous limit set by MINOS if v e appearance is ne- 
glected 12j, achieved from increased statistics and the 
reduction of systematic uncertainties in this analysis. 



tion algorithms allowing a reduction of the uncertainty 
due to the relative normalization between the two de- 
tectors from 4% to 2.2% 181. The absolute scale of 



Model 


#13 


x 2 /d.o.f. 


#23 


#24 


#34 


/- 


7714 = 7771 



11.5 


130.4/123 
128.5/123 


45.0± 7 
45.6±7 




0+ 17 
u - u -o.o 

n+ 25 

u - u -o.o 


0.22 
0.40 


7714 > 7713 




11.5 


130.4/122 
128.5/122 


45.0±7 
45.6±7 




n+ 17 

u - u -o.o 
n n +25 

u - u -0.0 


0.22 
0.40 



TABLE II: Best-fit values and uncertainty ranges in degrees 
for the angles of the two neutrino oscillation scenarios in- 
cluding a sterile neutrino. The results shown assume either 
no v e appearance or v e appearance at the MINOS 90% C.L. 
limit. The quantity f s is the maximum allowed fraction of 
disappearing that may transition to v s . 

The dominant systematic uncertainties were reevalu- 
ated, in some cases with more sensitive methodologies, 
leading to reductions from the previous analysis. Visual 
scanning techniques were used to asses the reconstruc- 



tive hadronic energy contributes a maximum of 10% to 
the uncertainty for E TCCO < 0.5 GeV to a minimum of 
6.5% when I£ r eco > 10 GeV. This new treatment com- 
bines a constant 5.6% uncertainty in the detector re- 
sponse to single hadrons with an energy-dependent un- 
certainty due to hadronization model and intranuclear 
effects [25j . The uncertainty on the relative energy scale 
between detectors, computed from interdetector calibra- 
tion studies, is 2.1%. The uncertainty due to ND-specific 
selections, determined by varying each selection crite- 
rion to assess its effects in the ND energy spectrum, 
is 10% for E rcco < lGeV, between 4.8% and 2.1% for 
0.5 < -E ra co < 2.5 GeV, and negligible for higher ener- 
gies The uncertainties stemming from FD-specific 
selections, evaluated in the same manner as the ND- 
specific ones, are calculated to be 5% for E ieco < 0.5 GeV, 
and between 2.5% and 1% for 0.5 < E ICCO < 120 GeV. 
The uncertainty in the size of the z^-CC background, 
computed usin g N D data acquired in different beam con- 
figurations [l2, 18 1, remains 15%. 

The effects of other systematic uncertainties, such as 
those due to the physics of neutrino interactions, largely 
cancel in the FD prediction due to the similarity in ma- 
terials and planar separation of the ND and FD. These 
additional uncertainties are the principal contributors to 
the error band shown in the ND energy spectrum in 
Fig. [D but their effects on the FD error band shown in 
Fig. [2] arc negligible. Table [ilTI summarizes the variation 
in best-fit values obtained for the oscillation parameters 
when the effects of the dominant systematic uncertain- 
ties are applied to the ND and FD energy spectra but not 
included in the fit. The largest systematic- induced shift 
in the value of R over the full energy range is 2.9% and 
is due to the relative normalization uncertainty. With 
the exception of the energy scale uncertainties, which do 
not modify the total number of events and thus have no 
effect on R, the other dominant systematic uncertainties 
each induce shifts of 2.3%. 

In summary, new results are presented from a search 
for active to sterile neutrino mixing, based upon a 
data sample with double the event statistics of previ- 
ous MINOS analyses. A total of 802 NC event candi- 
dates is observed in the FD data, compared to 754 ± 
28(stat) ± 37(syst) events expected from standard os- 
cillations. Therefore, no evidence for depletion of NC 
events is observed in the FD at a distance of 735 km from 
the production target. The most stringent constraint to 
date is placed on the fraction of active neutrinos that 
transition to sterile neutrinos, f s < 0.22 (0.40) at the 
90% C.L., where the number in parentheses denotes the 
limit assuming v e appearance. The results support the 
hypothesis that disappearance observed in MINOS is 
dominated by oscillations among active neutrino species. 
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TABLE III: Magnitude in degrees of the mixing-angle devia- 
tions introduced by the major systematic uncertainties from 
best-fit results in which systematic shifts have been neglected. 
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